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Abstract 

In the last approximately ten years, a number of codes using different approaches have been successful in solving the old 
problem of calculating the band gap of semiconductors with reasonable accuracy. The GW approximation (GWA), based on 
perturbation theory, some hybrid functional proposals, and the Wien2k code implemented with the modified Becke-Jonhson 
potential (mBJLDA) are compared in this work. In particular, we studied 27 semiconductors for which we found low 
temperature experimental data. We found that the GWA gives the most accurate predictions. The Wien2k (mBJLDA) code 
is slightly less precise, in general. The Hybrid functionals are less accurate, on the overall. The GWA is definitely the most 
precise existing method nowadays. In 88% of the semiconductors considered the error was less than 10%. Both, the GWA 
and the mBJLDA potential, reproduce the band gap of 15 of the 27 semiconductors considered with a 5% error or less. An 
(3JQ[ extra factor to be taken into account is the computational cost. If one would seek for precision without taking this factor into 
^ . account, the GWA is the method to use. If one would prefer to sacrifice a little the precision obtained against the savings 
in computational cost, the empirical mBJLDA potential seems the appropriate method. We include a graph that compares 
directly the performance of the best three methods, according to our analysis, for each of the 27 semiconductors studied. The 
situation is encouraging but the problem is not yet a closed issue. 

Short title: The band gap problem. 

PACS numbers: 71.15.Mb ; 71 .20.Mq; 71. 20. Nr 
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INTRODUCTION 



Khon and Sham equations [l[ are central to the practical application of Density Functional Theory (DFT) . To solve 
them, an approximation to the exchange and correlation energy is required from which an exchange and correlation 
potential is derived. The way in which this term is approximated is crucial to the proper description of the band 
structure of solids. The Local Density Approximation (LDA) the Generalized Gradient Approximation (GGA) 0- 
[H[ and the meta-GGA 0, Q, among others, describe very well the electronic band structure of even complicated 
metallic systems. They fail, nevertheless to account for the band gap value of semiconducting systems, a short come 
known for several years now Q. Efforts to solve this problem were done since long ago. Approximations as the "scissor 
operator" 0] , the Local Spin Density Approximation, LSDA+UJ0 and methods based on the use of Green's functions 
and perturbation theory as the GW approximation, GWA I3-12J, were proposed. In the last ten years, these efforts 



gave rise to substantially improved results. Some of the new proposals include, the screened hybrid functional of 
Heyd, Scuseria and Ernzerhof (HSE) (l3l4l5l] and the middle-range exchange and correlation hybrid functional of 
Henderson, Izmaylov, Scuseria and Savin (HISS) (lil . ■ Another recent proposal is the empirical potential the 
modified Becke- Jonhson potential (mB JLDA) proposed by Tran and Blaha |18| . This potential was introduced to the 
Wien2k code in 2010. 

THE MBJLDA POTENTIAL 

Recently, we made a detailed analysis of the mBJLDA potential based on the calculation of the electronic band 
structure of 41 semiconductors (20| . We found an important improvement in the predictions of the band gap as 
compared to experiment. The mBJLDA potential (l8| is a empirical potential of the form 
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where p a (r) is the spin dependent density of states. t a (r) is the kinetic energy density of the particles with spin a. 
and V x B J(r) is the Becke- Roussel potential (BR) (21(. The c stands for 
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a and j3 are free parameters. The Wien2k 2011 code defines a = —0.012 and j3 = 1.023 Bohr 1 / 2 . These values 
are general but certainly fixed experimenting with several cases. A particular feature of this potential is that a 
corresponding exchange and correlation energy term , E xc [p], such that the mBJLDA potential is obtained in the 
usual way, namely, V xc = SE xc [p]/Sp, is not possible. As a consequence, a consistent optimization procedure to obtain 
the lattice parameters, the Bulk modulus and its derivative with respect to pressure are not actually possible. This is a 
consequence of the empirical character of this potential. For that reason, Tran and Blaha have proposed the empirical 
alternative that prior to a band structure calculation with the mBJLDA potential, the lattice parameter is found from 
either a LDA or a GGA optimization procedure and the result introduced into the code to perform the band structure 
calculation of the semiconductor system. Such a procedure gives rise to quite improved results as compared to the 
previous version of the Wicn2k code. It is known that the LDA underestimates as a rule, the lattice parameters and, 
on the contrary, GGA overestimates them. We have explored the possibility of using the averaged value as the lattice 
parameter, a,A vg , where a^vg = (o-lda + ogga)/2. Here a^oA^GGA) is the lattice parameter obtained from an LDA 
(GGA) optimization procedure. When aAvg is used as input into the Wien2k code implemented with the mBJLDA 
potential, a better agreement of the band gap value with experiment is obtained as compared to the results with either 
a>LDA or cigga- So this procedure turns out to give better results than the one recommended by Tran and Blaha and 
its extra computational cost is relatively low. A surprising result was, nevertheless, obtained when the experimental 
low temperature lattice parameter, olt- was introduced instead. Unexpected deviations of the band gap value from 
experiment as big as 48% were obtained (20[. This is a disturbing result since the lattice parameters obtained from 
any optimization procedure are judged to be as good as the deviation from the experimental lattice parameter value 
is small, and so one expects to get the best result (the minimum deviation of the predicted band gap value from 
experiment) when the experimental lattice parameter is used. This is not the case. This fact throws doubts on the 
meaning of the optimization procedure altogether when the empirical mBJLDA code is employed. Nevertheless, we 
stress that the results obtained for the band gap value of semiconductors using the mBJLDA potential represents a 
relevant improvement at relatively low computational cost, a fact that we will emphasize below. 
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THE HSE METHOD 



Hybrid functionals are a linear combination of Hartree-Fock, LDA and GGA terms and were proposed initially 



with the aim of improving LDA and GGA in the calculation of the energy bands of molecules [22|, |23| . More recently, 
hybrid functionals were used as an effort to improve the old-standing problem of the band gap of semiconductors; 
they include the Heyd-Scuseria-Ernzerhof (HSE) functional [l3| proposed in 2003. It combines a screened short-range 
HF term and a screened short- and long-range functional proposed by Perdew, Burke and Ernzcrhof (PBE) The 
screened terms in HSE result from splitting the Coulomb operator into short- and long-range terms in the following 
way 

1 _ erfejur) erfjujr) 
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where the complementary error function erfc(ujr) = 1 — erf(u>r) and u> determines the range. The functional form of 
HSE is based on the hybrid functional of Perdew, Burke and Ernzerhof (PBEh) [24| (also known in the literature as 
PBE1PBE and PBEO) [HHU, as follows 

E PBEh = aE BP + (1 - a)E PBE + E PBE + E PBE , (4) 

The expression for the HSE exchange-correlation energy, E XC SE , is 

Eg SE = aE BPSR (u) + (1 - a)E^ PBE ' SR (^) + E^ PBE - LR (u) + E PBE , (5) 

where E BFSR is the HF short-range functional (SR), E% PBE > SR and E^ PBE ' LR are the short-range (SR) and long- 



range (LR) components of the PBE functional, a is a mixing constant that is derived from perturbation theory [27 1. 
In the literature, the functional HSE appears as HSE03 and HSE06. The difference is in the choice of the value of 
uj. We will refer to the HSE03 simply as HSE, in this work. In 2005, Heyd et al. 14| reported a study of the band 



gap and lattice parameters of semiconductor compounds using the HSE functional. We will comment on these results 
below. 

Recently, Marques et al. [28[ have proposed to relate the mixing constant a to dielectric properties of the solid. 
They took a = l/e p J 3E . Their calculation using the hybrid functional PBEO improves the predictions for the band 
gap value of the 21 semiconductors considered as compared to the original formulation. Furthermore, they used a ~ g, 
ya~j 4 where g takes the form of the term in parenthesis in eq.(2). They introduced this form of the a parameter 
into the hybrid functionals PBEO and HSE06, respectively, and got an improved result. These proposals improve the 
performance of the hybrid functionals at no extra cost. We will comment further on these results below. 



THE HISS POTENTIAL 

Another successful potential to calculate the band structure of semiconductors is the middle-range hybrid exchange 



and correlation Henderdon-Izmaylov-Scuderia-Savia functional (HISS) [lall2|- ^ a ^ so uses the PBE potential but in 
a different way, 

pHISS _ j^PBE , „ (jpMR-HF p MR-PBE\ ( R \ 
Arc -Arc +C M R{E X - E x ) (b) 

where the last two terms in parentheses are the middle-range (MR) exact exchange and middle-range PBE exchange 
energies, given by 

1 _ erfc{ui S Rr) + erf(u S Rr) + erf(uj S Rr) ~ erf(uj LR r) 
r ^ r / r s r / 

SR LR MR 

In 2012, Lucero et al. [29| reported their study of the band gap and lattice parameters of some semiconductor 
compounds using HISS, with oj,sr = 0.84a Q " 1 , lolr = 0.20<2q 1 y cmr = 0.60. These values were determined by fitting 
them to some atomization energies, barriers hights and values of the gap for some compounds. 
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THE GW APPROXIMATION (GWA) 



As it is well known, the many-body Shrodinger equation contains the Coulomb interaction term which is a two- 
body potential and creates the difficulty to solve it for realistic systems. To address this problem, the Hartree-Fock 
Approximation (HFA) adds to the average Coulomb potential (the Hartree term) a non-local exchange potential 
which reflects the Pauli Exclusion principle. The energy gap of semiconductors predicted in this way turns out to 
be in most cases too large. This is due to the neglect of correlations or screening which are crucial in solids. To 
simulate the effect of correlations, Slater introduced the Xa approximation which may be regarded as a precursor 
of the modern DFT. In DFT, the ground state energy can be proved to be a functional of the ground state density 
but the explicit form of the functional is not known. The minimization of the total energy functional with respect to 
the density gives the Kohn-Sham equations. The unknown exchange and correlation potential is approximated either 
by the local density approximation (LDA), the generalized gradient approximation (GGA) or the meta-GGA, among 
others, which describe metals well but fail to account for the band gap of semiconductors. The empirical mBJLDA 
potential is a response. An alternative way to deal with this problem is the GW approximation (GWA). It is derived 
from many-body perturbation theory [3(| ■ The form of the self-energy in the GWA is the same as in the HFA but the 
Coulomb interaction is dynamically screened remedying the most serious deficiency of the HFA. The corresponding 
self-energy is therefore non-local and energy dependent. The Green function is obtained from a Dyson equation of the 
form G — Gq + GqT,G where Go describes the direct propagation without the exchange and correlation interaction 
and £ contains all possible exchange and correlation interactions with the system that an electron can have in its 
propagation. The GWA may be regarded as a generalization of the HFA but with a dynamically screened Coulomb 
interaction. The non-local HFA is given by 

(occ) 

^(r.r'j^fW^rXr-r') (8) 

Where v(r — r') is the bare Coulomb interaction. The GWA corresponds to replacing the bare Coulomb interaction 
v by a screened interaction W. In the language of perturbation theory this corresponds to 

E x (r,r',w) = ^- [ du)'G{r,r',u+u>')W(r,r',u)') (9) 
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For details see ref. [Ill ] . We will analyze the results obtained for the semiconductor gap value using this approximation 
in what follows. 

In 2005, Rinke et al. [3l[ using the so called OEPx(cLDA)+GW approximation obtained a reasonable agreement 
with experiment when calculating the band gap of a certain number of semiconductors. In 2007, Shishkin et al. [32 



used a self-consistent GWA (GWA + DFT), and the self-consistent GW approximation with attractive electron- hole 
interaction, scGW(e-h) accounted quite well for the experimental band gap of several semiconductors. 

Now we proceed to analyze some of the different offers in the literature in what the calculation of the band gap of 
semiconductors is concerned compare their results among themselves. 

ANALYSIS 

In Table HI we present the results obtained for the band gap using different methods and approximations. The first 
two columns refer to the resulting values using LDA and mBJLDA [20| as in version 2011 of Wien2k code. We have 
used as lattice parameter the average of the values obtained from an LDA and a GGA optimization which, as we found 
in a previous work (20l |. gives the best results for the gap as compared to experiment when the mBJLDA potential is 
used. Next we report the values obtained with the hybrid HISS and HSE06 functional HSE (HQ and GWA. In 
the column denoted as GWA, we include the most precise predictions for the gap reported in the literature usin g ei ther 
the self consistent GW (scGWA) [33l - l40j or the scGWA with attractive electron- hole interaction, scGW(e-h) j32j ]. 
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TABLE I. We compare the results for the gap (E g ), in eV, obtained with the Wien2k(LDA) code, with mBJLDA potential, with 
the hybrid functionals HISS, HSE06, HSE and with the GWA (see text). The crystal structure and percentage difference with 
respect to the experiment is shown in parenthesis. The minus sign means that the calculation underestimates the experimental 
value. The experimental data are from refs. [41H49| |. 
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5.6 a (2.2%) 
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0.00 ( 
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5.00 
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6.50(-16.3%) 
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;-4i%) 
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;-46%) 
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2.45 


GaN(B3) 
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2.23 (57.0%) 


1.61 (13.4%) 
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1.44 d (1.4%) 


1.42 


AlSb(B3) 
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2.05 (22.0%) 
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1.64 d (-2.4%) 
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GaSb(B3) 
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[-91%) 


0.90 (9.8%) 


1.31 (59.8%) 
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(-41%) 
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(-69%) 
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0.38 
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3.81 
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1.19 


(-58%) 


2.74 (-2.8%) 


2.93 (3.9%) 


2.27 (-19.5%) 


2.32 (-17.7%) 


2.84 f (0.7%) 


2.82 
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1.20 


(-50%) 


2.38 (-0.4%) 


2.77 (15.9%) 


2.16 (-9.6%) 


2.19 (-8.4%) 


2.57 f (7.5%) 


2.39 


MgS(B3) 






5.18 (-4.1%)* 


5.17 (-4.3%) 


4.48 (-17.0%) 


4.78 (-11.5%) 




5.40 


MgTe(B3) 






3.59 (-2.2%)* 


3.91 (6.5 %) 


3.49 (-4.9%) 


3.74 (1.9%) 




3.67 


GaN(B4) 


2.06 


(-41%) 


3.13 (-10.6%) 


4.23 (20.9%) 


3.14 (-10.3%) 


3.21 (-8.3%) 


3.5 g (0.0%) 


3.50 


InN(B4) 


0.03 


(-96%) 


0.82 (15.5%) 


1.51 (112.7%) 


0.66 (-7.0%) 


0.71 (0.0%) 




0.71 


A1N(B4) 


4.11 


(-34%) 


5.53 (-10.7%) 


6.62 (6.9%) 


5.50 (-11.2%) 


6.45 (4.2%) 


5.8 s (-6.3%) 


6.19 


ZnO(B4) 


0.76 


(-78%) 


2.76 (-19.8%) 








3.2* (-7.0%) 


3.44 


A(%) 


60 


.2% 


8.4% 


34.1% 


11.5% 


10.2% 


5.7% 




The experimental band gap at room temperature 


BP(B3) 


1.15 


(-43%) 


1.83 (-8.5%) 


2.43 (21.5%) 


2.21 (10.5%) 


2.16 (8.0%) 




2.00 


BN(B3) 


4.39 


(-29%) 


5.85 (-5.6%) 


6.69 (7.90%) 


5.90 (-4.8%) 


5.99 (-3.4%) 


7.14 (15.2%) 


6.20 


MgSe(Bl) 


1.71 


(-31%) 


2.89 (17.0%) 


3.05 (23.5%) 


2.58 (4.5%) 


2.62 (6.1%) 




2.47 


BaS(Bl) 


1.93 


(-50%) 


3.31 (-14.7%) 


3.61 (-7.0%) 


3.21 (-17.3%) 


3.28 (-15.5%) 


3.92 (1.0%) h 


3.88 


BaSe(Bl) 


1.74 


(-51%) 


2.87 (-19.8%) 


3.14 (-12.3%) 


2.80 (-21.8%) 


2.87 (-19.8% 




3.58 


BaTe(Bl) 


1.37 


(-56%) 


2.24 (-27.3%) 


2.48 (-19.5%) 


2.22 (-27.9%) 


2.50 (-18.8%) 




3.08 


BAs(B3) 


1.23 


(-16%) 


1.72 (17.8%) 


2.14 (46.6%) 


1.89 (29.5%) 


1.92 (31.5%) 




1.46 


*with mBJ(flLD.4) 


t scGW(h-e) in Ref. [32]. 


t scGW in Ref. [33]. a Ref. [34]. 


b Ref. [35] . b Ref. 


[36]. 





d Ref. [37]. e Ref. [38]. f Ref. [39]. s Ref. [40]. h Ref. [50]. 

A(%) is the average of the absolute percent deviations. 

We calculate the percent deviation as follows Error(%) = (Eg eo ' — Eg Xpt ) 
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FIG. 1. The horizontal axis is the absolute percent error (\Error(%)\). Each horizontal line is drawn at the low-temperature 
experimental band gap value for each semiconductor. It compares, the three best results (GWA, mBJLDA, and HSE) according 
to our analysis. The inset shows in more detail the semiconductors with a gap between 2.2 to 2.5 eV for clarity.(see tableland 
text). 



In Fig. [U the band gap value is given in the vertical axis. Each horizontal line is drawn at the experimental low 
temperature band gap value for each of the 27 semiconductors considered. The horizontal axis represents the absolute 
percent error, (\Error(%)\) calculated as shown at the bottom of Table I. The three data on each horizontal line 
correspond to the result obtained using the mBJLDA potential, the HSE method and the GWA. So, for a particular 
semiconductor, the graph compares directly the performance of each of the three best methods as found in this work, 
(see table [XJ) . As it is very well known Q and as it appears in Table |H the LDA does not reproduce the experimental 
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values of the band gap of semiconductors. Furthermore, The Wien2k (LDA) code produces for MgS and MgTe a band 
structure which shows a direct band gap in contradiction with experiment pol ] . 

The results of the predictions obtained with the GWA are the most accurate with an averaged error of 5.7%. The 
empirical mBJLDA potential produces results with an averaged error of 8.4%. Next, the errors obtained with the HSE 
potential result in an averaged error of 10.2%, HSE06 (11.5%) and HISS (34.1%). The GWA, the mBJLDA potential, 
and HSE functional do better than the ones reported by Marques et al. [28| . They get results with averaged errors 
16.5%, 14.4% y 10.4% using the hybrid functional PBEQ too , PBEQ mix and HSE06 mix , respectively. In this paper, 
the authors suggest the interesting possibility that the mixing parameter should be related to physical variables. The 
performance of the GWA is highly accurate, 88% of the calculated results recorded here show less than 10% error. 
This is to be compared to the one obtained when using mBJLDA (74%), HSE (54%), HSE (42%), HISS (35%). On 
the other hand, more than a 20% deviation from experiment occurs when using mBJLDA in 7% of the cases, and with 
the GWA in 8% of the semiconductors studied, which is to be compared with HSE (15%), HISS (46%). All together, 
the best results are obtained when using the GWA; mBJLDA is next, but also HSE gives results with acceptable 
accuracy. A special case is the very-low-gap InSb. In this case none of the methods give less than a 20% error although 
it would be more reasonable to judge these results from the absolute deviation in electron-volts rather than from the 
percent deviation (see Table I). It is important for the overall picture to stress that the GWA and calculations with 
the mBJLDA potential present deviations less than 5% in 15 of the 27 semiconductors considered. When the HSE 
potential is considered, 12 of the 27 present less than a 5% deviation . 

One more observation. In the previous analysis we took into account only low temperature band gap data. In Table 
I, we also present some calculations for which we did not found experimental reports at low temperature. Since the 
calculations are done at OK, room-temperature measurements require extrapolation either using Varshni's law or a 
quadratic fit or any other suitable method which, in any case, generates an extra incertitude in the obtained OK data. 
If we rather use the high temperature data, the HSE potentials give a better agreement with experiment. 



CONCLUSIONS 

The accurate calculation of the band gap of semiconductors is a difficult task that has been the object of intense 
research with the result of important progress during the last approximately ten years. We presented in this work 
our analysis of some different solutions and compare their results among them. A group of 27 semiconductors (see 
Table I) for which we found low temperature data on the band gap value were considered. The results of the GWA, 
the Wien2k implemented with the mBJLDA potential and codes using a hybrid functional, HSE, and HISS were 
taken into consideration. The results reported by Marques et al. [28| were found to be less accurate than the ones 
of the GWA, the mBJLDA potential and the HSE functional (see text above for the precise definition used here). 
The GWA was found to give, all together the best results. The mBJLDA potential produces results slightly less 
accurate and HSE comes next. The two first methods give quite good results (prediction better than 5% for 15 
of the 27 semiconductors studied). In Fig. [TJ we compare the performance of the three best methods found in this 
analysis for each of the 27 semiconductors separately. It is important to stress the empirical character of the mBJLDA 
potential because it prevents the consistent definition of the optimization procedure which contrasts with the sound 
bases of the GWA. Even with the several theoretical non-properly solved issues, the mBJLDA potential gives rise to 
acceptable predictions of the band gap value as compared to experiment. An extra factor to be taken into account 
is the computational cost. If one would seek for precision without taking this factor into account, the GWA is the 
method to use. If one would prefer to sacrifice a little the precision obtained against the savings in computational 
cost, the mBJLDA potential seems the appropriate method. In conclusion, we can typify the state of matters with 
respect to the calculation of the band gap of semiconductors as follows. A quite precise method does exist, the 
GWA approximation. It's computational cost is higher. A relatively quicker code, the Wicn2k implemented with 
the mBJLDA potential, gives somehow less accurate results but quite acceptable at lower computational cost. Other 
methods do exist but are less accurate. Very recently [5~0 ]. the new approximation announced in the ref [EH was 
implemented in the Wien2k 12.1 code for public use. The new hybrid functional YS-PBE0 is "equivalent" to the HSE 
one, according to the authors. We will study this new functional in future work. 
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